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We establish a correspondence between the electric dipole matrix elements of a polyatomic symmetric top 
molecule in a precessing state and the magnetic dipole matrix elements of a magnetic dipole associated with 
an elemental spin F. It is shown that this correspondence makes it possible to perform quantum simulation 
of the single-particle spectrum and the dipole-dipole interactions of magnetic dipoles in a static external 
magnetic field B with symmetric top molecules subject to a static external electric field Edc- We further 
show that no such correspondence exists for *E molecules in static fields, such as the alkali metal dimers. 
The effective spin angular momentum of the simulated magnetic dipole corresponds to the rotational angu- 
lar momentum of the symmetric top molecule, and so quantum simulation of arbitrarily large integer spins 
is possible. Further, taking the molecule CH3F as an example, we show that the characteristic dipole-dipole 
interaction energies of the simulated magnetic dipole are a factor of 620, 600, and 310 larger than for the 
highly magnetic atoms Chromium, Erbium, and Dysprosium, respectively. We present several applications 
of our correspondence for many-body physics, including long-range and anisotropic spin models with ar- 
bitrary integer spin S using symmetric top molecules in optical lattices, quantum simulation of molecular 
magnets, and spontaneous demagnetization of Bose-Einstein condensates due to dipole-dipole interactions. 
Our results are expected to be relevant as cold symmetric top molecules reach quantum degeneracy through 
Stark deceleration and opto-electrical cooling. 



1 Introduction 

The realization that the phenomena of electricity and magnetism are related manifestations of the same 
field is a landmark in the history of physics, marking the first successful unification of theories. In spite 
of the fact that electricity and magnetism are derived from the same physical source, electric phenomena 
and magnetic phenomena are not identical. This can be seen classically by the fact that the electric field is 
a vector, transforming into its mirror image under inversion of all coordinates, while the magnetic field is 
a pseudo-vector, transforming into minus its mirror image under inversion of all coordinates. In quantum 
mechanics, this difference between electric and magnetic fields has the consequence that matrix elements 
of the electric dipole operator vanish unless the two states have opposite parity, while magnetic dipole 
matrix elements connect states which have the same parity. Hence, an elemental object can possess a 
magnetic dipole but not an electric dipole, as possession of an electric dipole would violate both parity and 
time -reversal symmetrjQ 

Even in light of the fundamental differences between electric and magnetic dipoles, interactions between 
two electric or magnetic dipoles are described by a Hamiltonian of the same form, namely |4| 
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1 As discussed also in Sec. [5] parity and time-reversal are not exact symmetries of the universe within the understanding of the 
standard model. However, their violations are very small. Ultracold molecules, incidentally, provide a highly accurate platform to set 
bounds on these violations fTH2]|3). 
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where T>i is the dipole moment operator T> of the i th particle with T> = d, fi for electric and magnetic 
dipoles, respectively. Also, r is the distance between the two dipoles, e r is a unit vector pointing from 
dipole 1 to dipole 2, and the coupling coefficients for the electric and magnetic dipole-dipole interactions 
are given by 

C d = 1/eo, = mo, (2) 

with e the vacuum permittivity and /i the vacuum permeability. The dipole-dipole interaction Eq. ([Tji is 
of interest from the point of view of many-body physics because it is long-range and anisotropic. The dif- 
ferences in coupling coefficients, Eq. Q, demonstrate a practical difference between electric and magnetic 
dipoles. Namely, the interactions between electric dipoles on typical atomic scales are roughly 10 4 times 
larger than for magnetic dipoles. Hence, quantum simulation of a gas of magnetic dipoles with electric 
dipoles would provide access to the physics of magnetic dipoles in more strongly correlated regimes, on 
shorter timescales, and at lower density than for the true magnetic dipoles encountered in nature. In this 
work, we show that such a quantum simulation is indeed possible. 

It is only in the last ten years that quantum degenerate gases with significant dipole-dipole interac- 
tions have been an experimental reality. The first strongly dipolar gas to reach quantum degeneracy was 
Chromium |3), which has a magnetic dipole moment of 6/is, where /is is the Bohr magneton. A clear 
manifestation of the effects of dipole-dipole interactions were seen in the <i-wave symmetry of the atomic 
cloud expansion following collapse of a Chromium Bose-Einstein condensate [6]. Even for the alkali 
metals with magnetic dipole moments of dipolar effects have been seen in some cases 0O. For 
Chromium, the ratio of the dipole-dipole interaction energy to the interaction energy due to isotropic col- 
lisions is 0.15, and so the effect of dipole-dipole interactions in this system is considered to be a small 
perturbation. Recently, degenerate gases of Erbium [9 1 and both bosonic [ 10] and fermionic [ 1 1 1 Dyspro- 
sium have been produced, which have magnetic dipole moments of 7/ib and 10/is, respectively. Assuming 
that the scattering lengths of Erbium and Dysprosium are both roughly 100 Bohr radii, as is the case for 
Chromium, this gives ratios of the dipole-dipole interaction energy to the isotropic interaction energy of 
0.67 and 1.33, respectively. Significant progress has also been made on the production of ultracold polar 
molecules [ 12], with KRb IPT31 being the first to reach the ultracold regime. Several other species are cur- 
rently being investigated [ 14 1, the majority of which are diatomic molecules consisting of two alkali metal 
atoms lfl"5l [TBI [T7l [T8l . The alkali metal dimers are 1 S molecules, having no electronic spin or orbital 
angular momentum, and so their angular momentum structure corresponds to that in Fig.[JJa). Namely, the 
total angular momentum J consists only of rotation, and is perpendicular to the internuclear axis. 

In many theoretical works on the many-body physics of dipoles, it is presumed that only the magnitude 
of the dipole moment has importance. In particular, it is supposed that the origin of the dipole moment, 
whether it is magnetic or electric, is of no importance. For the case of polarized dipoles, in which only a 
single quantum state with a resonant dipole moment is populated, this presumption is indeed true, as we 



discuss in Sec. 4.1 However, for the multi-component case, the structure of the dipole matrix elements 
within the allowed state manifold and the energies of these states are often very different depending on 
the origin of the dipole moments. In particular, while the alkali metal dimers typically have characteristic 
dipole-dipole interaction energies several orders of magnitude larger than for atoms with magnetic dipole 
moments, the energies and expected dipole moments of these molecules are very different from magnetic 
atoms. Magnetic dipoles associated with an elemental angular momentum F display a linear coupling to 
an externally applied magnetic field such that differences in energies between states with projections of F 
along the field direction of M p and Mj? ± 1 are the same. In addition, the magnetic dipole moment behaves 
as though it points along F, and so the matrix elements of the magnetic dipole moment along any direction 
in space can be non-vanishing. Loosely speaking, such a magnetic dipole moment is "always on." On 
the other hand, the electric dipole moment of a X S points along the inter-nuclear axis, and its expectation 
vanishes along any direction in space due to effective averaging by the rotation of the molecule. A large 
static electric field must be applied in order to polarize the dipole moment along a particular direction 
in space. However, as the coupling of a 1 S molecule to an external electric field occurs only in second 
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Fig. 1 Schematic of angular momentum geometry for (a) KRb, a linear rigid rotor; (b) CH3F, a polyatomic symmetric 
top molecule; and (c) Cr, a highly magnetic atom. Here, we display both the space-fixed Z-axis and the molecule-fixed 
2-axis as well as the projections of the angular momentum along them, see the discussion following Eq. d3). In the 
molecular cases, panels (a) and (b), the total angular momentum consists only of rotation. For the linear rigid rotor 
in panel (a) this angular momentum has no projection on the molecular symmetry axis, but for the symmetric top 
molecule in panel (b) the angular momentum can have some projection K on this axis. 



and higher orders of perturbation theory, the energy levels and dipole moments do not resemble those of 
a magnetic dipole. Thus experiments like those of Bruno Laburthe-Tolra on the many body physics of 
Chromium atoms in optical lattices lfT9l l20l |2T1 l22l cannot be reproduced or explored in the many 1 E 
molecules now under intensive development in laboratories all over the world ifPfl [T31 [T6l [171 [181 . 

We will overcome these difficulties in quantum simulation of magnetic dipoles with 1 S molecules by 
focusing on symmetric top molecules. A symmetric top molecule is a polyatomic molecule with cylindrical 
symmetry, resulting in a doubly degenerate eigenvalue of the inertia tensor. This implies that the angular 
momentum of a symmetric top molecule can have some component in the direction of the symmetry 
axis, see Fig. [TJb). In particular, when the projection of the angular momentum along the molecular 
symmetry axis is non-vanishing, a symmetric top molecule displays a linear response to an externally 
applied electric field, so that the energies and dipole matrix elements mimic those of an elemental magnetic 
dipole. In addition to their fundamental interest in quantum simulation given in this paper, symmetric 
top molecules are also promising candidates for reaching quantum degeneracy precisely because of their 
sensitivity to electric fields. This sensitivity allows for molecules to be slowed by the application of an 
electric field gradient known as a Stark decelerator 11231 . In addition, opto-electrical cooling, where electric 
field interaction energy is used rather than photon recoil to remove energy from translational motion, is now 
a demonstrated technology for the symmetric top molecule CH 3 F [24). For this reason we take CH 3 F as 
the principle physical example of a symmetric top molecule for this article. To date, the most successful 
means of producing ultracold molecules have been assembly of ultracold molecules from ultracold atoms 
by magnetic lfT3ll or optical [25] means or direct laser cooling of molecules [26|, all of which work only 
for specific species. In contrast, opto-electrical cooling has the potential to work for any symmetric top 
molecule, of which CH 3 F represents an initial experimental demonstration. Opto-electrical cooling hence 
opens up a new paradigm for the production of ultracold quantum gases, especially of common molecules 
relevant to chemistry. 

The ability to mimic the single-particle and two-particle matrix elements of quantum magnets with sym- 
metric top molecules also implies the ability to engineer quantum simulations of dilute gases or crystals of 
quantum magnets with symmetric top molecules. These quantum simulations open up a wealth of many- 
body physics, presented in Sec.[4j including new regimes of many-body physics which are not accessible 
with naturally occurring magnetic dipoles. For example, we show that polarized gases of symmetric top 
molecules in a single quantum state interact much more strongly than single-component gases of magnetic 
dipoles, and require much smaller electric fields than 1 E molecules. Long-range and anisotropic lattice 
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spin models can be simulated with symmetric top molecules in deep optical lattices. In addition to the 
long-range version of the XXZ model where simulation is also possible with linear rigid rotors 11271 l28ll . 
we demonstrate how to simulate spin models which do not conserve magnetization, a feat which requires 
fine-tuning of a large number of fields to achieve with linear rigid rotors |29 1. A major advantage of simu- 
lations with symmetric top molecules is that any effective spin may be chosen without fine-tuning of fields 
by populating a different rotational manifold. Furthermore, the exchange statistics of the molecule is not 
associated with this effective spin, and so symmetric top molecules provide a quantum simulation of a 
fermionic particle with a large integer spin, something which does not occur for actual magnetic dipoles. 
Using the connection of large-spin physics, we discuss using symmetric top molecules to simulate molec- 
ular magnets in condensed matter. Finally, we discuss simulations of the spontaneous demagnetization of 
degenerate Bose gases such as has been seen for Chromium ET1I221 . Here, simulations with symmetric top 
molecules not offer the advantages of shorter timescales for dynamics and lower requirements for density, 
but also allow the possibility of exploring new strongly correlated regimes which would not be accessible 
to any magnetic atom. 

This paper is outlined as follows. In Sec.|2]we review the rotational structure of symmetric top molecules 
and discuss the coupling of symmetric top molecules to external fields. In Sec. [3] we review the structure of 
magnetic dipoles and present a correspondence between the energies and dipole matrix elements of sym- 
metric top molecules and magnetic dipoles. Sec. |3.2| in particular Table [3] contains such correspondence, 
which is the main focus of this paper. In Sec.[4]we show this mapping can be used to explore many-body 
physics, and the advantage of using symmetric top molecules versus actual magnetic dipoles. In Sec. [5] 
we discuss the influence of hyperfine structure on our results, focusing on the case of CH 3 F as a typical 
symmetric top molecule. Finally, in Sec. [6] we conclude. There are three appendices. Appendix A contains 
some useful matrix elements for symmetric top molecules. Appendix B discusses the relationship between 
the microscopic angular momentum structure of a highly magnetic atom and the model of an elemental 
magnetic dipole used in the main text. Finally, Appendix C discusses the mapping between symmetric top 
molecules and magnetic dipoles presented in Sec. |3.2| in second quantization. We hope that this makes our 
results accessible not only to researchers in molecular physics, but also to those in condensed matter. 

2 Rotational structure of symmetric top molecules and coupling to exter- 
nal fields 

2.1 Quantum mechanics of symmetric top molecules 

For quantum gases at ultralow temperatures thermal energy is insufficient to create excitations in the struc- 
tural or electronic degrees of freedom of the atoms or molecules making up the gas. Hence, in the theory 
of ultracold atomic gases, modeling atoms as point-like particles with fixed angular momentum degrees 
of freedom is appropriate. For molecules in low-lying rotational states, we assume that molecules have 
the equilibrium spatial structures of their nuclei. Additionally, we can assume that the molecules remain 
always in their lowest vibrational and electronic states. Hence, the rotational motions of nuclei in closed- 
shell molecules are well modeled by the quantum mechanics of rigid bodies, known as the rigid rotor 
approximation (RRA) QUI . Within the RRA, the rotational energy of a molecule is given by the Hamilto- 
nian 

J2 J2 j2 

H ^ = 2t + 2t + 2t< (3) 

where x, y, and z are spatial directions along the principal axes of inertia of the rigid body in the body-fixed 
frame, or molecule-fixed frame which rotates with the molecule. We shall also make use of the space-fixed 
frame of coordinates X, Y, and Z, which is the usual system of laboratory coordinates. Also in Eq.([3]l, J x , 
J y , and J z are projections of the angular momentum J along the x, y, and z directions, respectively, and 
I x , I y , and I z are the principal moments of inertia. For a symmetric top molecule, the inertia tensor has a 
doubly degenerate eigenvalue. Defining the z-direction as lying along the symmetry axis of the molecule, 
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we have that I x = I y = Ij_. Hence, the rotational Hamiltonian for a symmetric top molecule can be 
written as 

4ot = B Q J 2 + (A - B ) J 2 Z , (4) 
where we have defined two rotational constants as 

The subscript "0" in the rotational constant sets Bq apart from the strength of the applied magnetic field 
B, and reminds us that this rotational constant refers to the ground vibrational level of the molecule. The 
square of the total angular momentum J 2 is quantized with eigenvalues J (J + 1), where the values J 
are non-negative integers, and the projection along a molecule-fixed quantization axis, which we take to 
be the z-direction, is also quantized with integer eigenvalue K. For a general symmetric top molecule, 
the projection K can take on any integer value in the range —J<K<J. In addition, the angular 
momentum is also quantized along the space-fixed Z-direction, see Fig. |TJb). The associated eigenvalue 
of Jz is denoted by M, and takes on integral values in — J < M < J. These three quantum numbers J, 
K, and M completely specify the eigenstates of the symmetric top rigid rotor Hamiltonian Eq. Q. The 
corresponding eigenfunctions are given explicitly in the Euler angle representation by 



(u\JKM) = \/ 2 ^^», (6) 
with corresponding eigenenergies, 

Ejkm = B J(J+l) + (A -B Q )K 2 , (7) 

where &f^ K (ui) are the matrix elements of the Wigner 7?-matrix [ 30 1 that transforms the space-fixed frame 
onto the molecule-fixed frame by three Euler angles oj = (cf>, 8, x). Throughout this paper, the appearance 
of w in any TJ-matrix or spherical harmonic refers to the three Euler angles rotating the molecule-fixed 
frame to the space-fixed frame and not to an angular frequency. States with the same \K\ ^ are doubly 
degenerate, due to the cylindrical symmetry of the Hamiltonian Eq. Q. For a discussion of the degree to 
which the ±\K\ degeneracy persists for symmetric top molecules beyond the RRA, see Sec. [5] 

For the case of a linear rigid rotor, the rotational angular momentum must be perpendicular to the 
symmetry axis of the molecule, and thus its projection K along the symmetry axis is identically zero, as 
shown Fig.[jja). In this case, the wave functions become proportional to spherical harmonics Yjm(9, 4>) 

(u\J0M) = *Y JM {6,<f>), (8) 

V 2,TT 

with corresponding eigenenergies, 

Ejom = B J(J + 1). (9) 

The disappearance of the angle \ from the wavefunction Eq. ([8| is a consequence of the linear structure. 
The two spectra Eq. |7]) and Eq. |9]) and their associated degeneracies are compared in Fig. [2] 

It might seem strange that we can specify two projections of an angular momentum, the projection 
along the space-fixed Z-axis M and the projection on the molecule-fixed z-axis K, as it is known that 
one cannot specify the projections of an angular momentum along two space-fixed axes simultaneously. 
The reason that we can specify these two quantum numbers, or, equivalently, the reason why Jz an d 
J z commute, is that the angular momentum operator J expressed in the molecule frame is not a true 
angular momentum operator in the sense of generating of a Lie group. Rather, the molecule-fixed angular 
momentum operator obeys anomalous commutation relations in which the Lie bracket has an extra minus 
sign, i.e., [J x ,J y ] — —iJ z . These anomalous commutation relations can be viewed as time-reversed 
rotation as seen from the frame which rotates with the molecule [31 1. 
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Fig. 2 Rotational energy spectra of (a) linear rigid rotors Ejom and (b) symmetric top molecules Ejkm- The 
degeneracy of a given level is denoted by the number to the upper right of the level. The dashed ellipses in (b) are 
guides to the eye of fixing J and varying \K\, called fixed J-manifolds. For simplicity, panel (b) assumes a prolate top 
with Aq > Bo- However, the scale of Ao/Bq for a given symmetric top molecule may be very different than the one 
displayed. 



2.2 The Stark effect and matrix elements of the dipole moment 

Applying a DC electric field to polar molecules causes the Stark effect, which is described perturbatively 
via an interaction term between the dipole moment operator of molecules d and applied DC electric field 
Erjc> that is, 

#Stark = — d • EdC • (10) 

In the following we assume that the DC electric field is applied along the space-fixed Z-direction with 
a strength of -Edcn i- e -> Edc = i?DC e z- F° r symmetric top molecules in the RRA, first-order perturbation 
theory yields the energy shift as 

AEjKM = {JKM \ -ffstark \ JKM) 

where d is the permanent electric dipole moment of the molecule, see Appendix A. For a fixed J, the 
(2 J + l)-fold degeneracy indexed by M has been lifted by the Stark effect, but there still remains a 2-fold 
degeneracy in the product of K and M when \K\ > 0. For a symmetric top molecule there is a linear 
Stark effect whenever \K\ ^ 0. That is, the Stark effect is linear both in Edc an d M for states with the 
same J and K, i.e., a fixed ( J, i\~)-manifold. We will call states with \K\ ^ precessing states, due to 
the fact that classically a state with if / precesses rather than tumbles ll32l . This precession is what 
gives a nonzero expectation of the dipole moment along the field direction and, hence, the first-order Stark 
effect. In addition, as the matrix elements Eq. ( fTT) i are the only ones which are non-vanishing for a given 
J and \K\, the states | J KM) are the eigenstates of the rotation and Stark Hamiltonian up to first order in 
dE-oc/Bo. For discussion of dipole matrix elements, it is convenient to define the space-fixed spherical 
basis {e p } p= o ) ±i, given by the unit vectors QUI 



e±i = T(ex ±iey)/V2, e = e z ■ (12) 

The expectation of the dipole operator along space-fixed spherical direction e p , d p = d • e p , in the basis 
\JKM) with J fixed, is 



(JK'M'\d p \JKM) = d(-l) J - M '( _ J M , l p m) K J J(j\\) 5kk '> (13) 
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where ( ;;; ) denotes the Wigner 3-j symbol [ 30 1 . The dipole matrix elements can be written in the form of 
the Wigner-Eckart theorem for a spherical tensor operator 



J 1 J 



(JK'M'\d p \JKM) = (-1) J - M I _ M , p " M \ (JK'\\d\\JK) , (14) 
via identification of the reduced matrix element 

(JK>\\d\JK) = dK^j Skk>. (15) 

There are three important implications of the linear Stark effect for our purposes. First, according to 



Eq. (Ill, the energetic separation between states | J KM ± 1) and | J KM) has the same magnitude for 
all values of M. Second, for any non-vanishing electric field strength -Edo we obtain an expected dipole 
moment along the space-fixed Z-axis of dKM/[J(J + 1)], which becomes half the permanent dipole 



moment when J = KM = 1. Finally, due to the spherical tensor structure in Eq. ( 14 1, the dipole operator 
behaves as though it simply points in the direction of J when our attention is restricted to a single (J, K)- 
manifold. This mimics the behavior of magnetic dipoles, as will be discussed in Sec. |3. 1| These features 
of symmetric top molecules are compared with the behavior of a linear rigid rotor in Fig. [3] 

We now contrast these results for symmetric top molecules with the results for a linear rigid rotor. For 



a linear rigid rotor the first-order perturbation Eq. (Hi identically vanishes for any rotational eigenstate 



since K = 0. Then, the second-order Stark shift is given as 

vp (2) \(JKM\H St ^\J'K'M>) | 2 

l) 2 -K 2 ] [(J+l) 2 -M 2 ] (J 2 -K 2 )(J 2 - 




+ 1) 3 (2J + l)(2J + 3) J 3 (2J + 1)(2J- 1 




which yields for K = 0, 



AF (2) (dE DC ) 2 f J(J+l)-MI 2 ) 

^jom - 2Bq \ j(j + i)( 2 J _l)(2 J + 3) / ■ ( ' 

Here, pairs of states with the same \M\ remain degenerate. The expected dipole moment of a linear rigid 
rotor along the Z-direction can be obtained perturbatively from the Feynman-Hellman theorem 11391 as 

As opposed to the space-fixed dipole moments of a symmetric top Eq. ( fl"3"j ), which are independent of -Edc> 
the dipole moments of a linear rigid rotor depend linearly on -Edcn an d so are small in fields which are 
small compared to Bq / d. On the other hand, for fields which are large compared to B^/d the expected 
dipole moments approach the permanent value d, while states with different values of \M\ are separated by 
energies of the order of the rotational constant Bq. Because a linear rigid rotor requires non-perturbative 
fields in order to access a significant fraction of the permanent dipole moment along the field direction, the 
eigenstates of the rotational and Stark Hamiltonian for rigid rotors are no longer well-approximated by the 
rotational eigenstates | JOM) . However, due to the facts that there are no crossings between states which 
correspond to different rotational quanta J in zero field and that M is a good quantum number even in the 
presence of Edc, we can label the states as | JOM) , where J is an adiabatic label correlating to zero field. 

The behavior of linear rigid rotors and symmetric top molecules in a static electric field are contrasted 
in Tables [T] and [2] and Fig. [3] The parameters in Table [T] are obtained from experiment for KRb, but are 
computed from density functional theory or other ab initio means for the other species. In contrast, the 
parameters in Table [2] are obtained from experiment due to the ready availability of these symmetric top 
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Table 1 The exchange statistics (f=fermion, b=boson), rotational constants Bo, permanent dipole moments d, and 
polarizing fields E p such that dE c / Bq — 6 for several relevant X E alkali metal dimers 1 33 . 34 35|. At the field E c , 
the expected dipole moment of the ground state along the field direction obtains roughly 70% of its permanent value. 
The large fields E c required to polarize the dipole moment are in contrast to symmetric tops, in which essentially no 
field is required to polarize a large fraction of the dipole moment. 





Statistics 


B (GHz) 


d (Deybe) 


E p (kV/cm) 


6 Li 133 Cs 


f 


6.520 


5.52 


14.07 


23 Na 40 K 


f 


2.826 


2.76 


12.20 


B/ Rb 133 C S 


b 


0.504 


1.25 


4.80 


4U K »Y Rb 


f 


1.114 


0.566 


23.4 


6 Li 23 Na 


f 


12.735 


0.56 


271.043 



Table 2 The exchange statistics, rotational constants Bo, permanent dipole moments d, and critical fields E c such 
that dE c /Bo — 1 for selected symmetric top molecules |36 37 38 1. The critical field is an estimate for significant 
breakdown of the linear Stark effect regime. Thus, a tremendous range of electric fields can be used in experiments 
without going beyond the linear Stark approximation. 





Statistics 


B (GHz) 


d (Deybe) 


E c (kV/cm) 


12 CH 3 F 


b 


25.536 


1.850 


27.419 


la CH 3 F 


f 


24.862 


1.850 


26.696 


CH 3 C1 


b 


13.292 


1.87 


14.120 


CH3I 


b 


7.501 


1.62 


9.19 


CH 3 CN 


f 


9.1988 


3.92 


4.661 



species. As shown in Table [T] the fields required to significantly orient the ground state of alkali dimer 
molecules along the field direction via the second-order Stark effect range from a few to a few hundred 
kV/cm. In contrast, due to the linear Stark shift, the symmetric top molecules given in Table|2]obtain half of 
their permanent dipole moment for the 1 1 1 — 1) state in any non-vanshing field. The critical fields E c given 
in Table [2] represent a rough estimate where the purely linear Stark effect, i. e. , first order perturbation 
theory, becomes of the same order of magnitude as the second-order Stark effect. A comparison between 
the energetic spectrum and dipole moments of a symmetric top with J = K = 1, a linear rigid rotor 
with J = 1, and a magnetic atom with F — 1 are given in Fig. [3] While panels (a) and (b) of Fig. [3] 
are numerical data from non-perturbative calculations, panels (c)-(f) are schematics of the behavior where 
linear coupling to the applied field holds, see Tables [T] and [2] Figure [3] stresses the similarity between 
symmetric top molecules and magnetic dipoles both in terms of energies and dipole matrix elements, and 
also stresses the difference between linear rigid rotors and magnetic dipoles. 

2.3 Interaction of symmetric top molecules with other external fields 

In this section we discuss two other means by which symmetric top molecules couple to external fields. 
The first is the coupling of symmetric top molecules to off-resonant optical potentials. Optical trapping 
is important for our proposal, as it is the only means by which we can trap both the high-field and low- 
field seeking states contained in a particular (J, i^)-manifold simultaneously. The second topic is the 
coupling of magnetic fields to the magnetic moment generated by rotation of the molecule, known as the 
rotational Zeeman effect. Magnetic fields are an important means of separating out hyperfine structure 
such that single quantum states may be addressed |34|, and so it is also important to understand the effects 
of magnetic fields on the rotational structure of the molecule. The Zeeman effect due to intrinsic nuclear 
magnetic moments and associated nuclear hyperfine structure will be discussed in Sec. [5] 
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d a (0,±l)__J 
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dEv C /B„ 6 




Electric field strength 



2 |d s (l,±l)| |dr(l)| 



v 
o 

s 



d«(l,0) 
Electric field strength 




Magnetic field strength 

^ IM±i)l \iit\ 



Magnetic field strength 



Fig. 3 Energies and dipole matrix elements of a linear rigid rotor in the states jlOAf) (panels (a)-(b)), a symmetric 
top molecule in the states jllA/) (panels (c)-(d)), and a magnetic atom in the states \F = 1, Mf) (panels (e)-(f)). The 
behavior of the energies and electric dipole matrix elements of a symmetric top molecule in an electric field match the 
behavior of the energies and magnetic dipole matrix elements of a magnetic dipole in a magnetic field. The behavior 
of a linear rigid rotor in an electric field is significantly different. The dipole matrix elements are given in terms of 
the transition matrix elements dr(K) = (IK, ±l|d±i|l-K"0) (/it = (1, ±l|A±i 1 10} for the magnetic dipole) and the 
resonant matrix elements d R (K,M) = (iKM\d \lKM) (v-r(M f ) = {lM F \fi \ IM F ) for the magnetic dipole). 
See Sec. |3.1| for a complete discussion of magnetic dipole matrix elements. 



Let us consider the application of an optical field 

E opt (r,t) = E opt (r) e"^* + E* pt (r) e^* , (19) 

where the frequency of the field Lu opt is far detuned from any molecular resonances. We will suppose that 
the field separates into spatial and polarization components as 

l 

E op t (r) = E opt (r) £ e p e* p , (20) 



where e p are the spherical basis vectors Eq. ( 12 1 and the polarization vector e has unit norm, e* ■ e = 1. 



When the field is far off -resonant, the coupling of the molecule to the field may be described as 

H opt = -E* pt (r) • a (w opt ) • E opt (r) , (21) 

where a(uj opt ) is the dynamical polarizability tensor ll40l evaluated at the field frequency Lu opt . The 
symmetry of the polarizability tensor is the same as the inertia tensor. That is, a (w op t) also has a doubly 
degenerate eigenvalue, and the principal axes of the two tensors are parallel. It should be noted, however, 
that the polarizability tensor need not be prolate even if the inertia tensor is; only the symmetry is specified 
a priori. The cylindrical symmetry implies that the polarizability tensor may be written in the molecule- 
fixed spherical basis {e g } 9= o,±i as 

a (w op t) = ckji (w pt) e 9=0 ® e g=0 + «_l (w op t) [e g=i ® e* =1 + e g= _i ® e* q= _ 1 ] , (22) 

where the molecule-fixed z-axis is defined to be along the symmetry axis. Transforming from the molecule- 
fixed frame indexed by q to the space-fixed frame indexed by p via the transformation BTIl 

l 

q=-l 
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with lj the Euler angles connecting the two frames, the Hamiltonian Eq. (21 1 may be written 

2 

H opt = -|£ op t(r)| 2 <iK pt ) + Aa(u; pt) C p ] ^) T v 

p=-2 

Here we have defined products of polarization components T p as 



T 



±2 



and unnormalized spherical harmonics C p as 



~e T ie* ±l , T ± i — — (e a e* ±1 - £=fi £ o) ■ T o - ko| 2 - 1/3 , 



-in 



■Y ep {9,<t>) , 



2^+1 

Also, we introduced the polarizability tensor invariants 

2a± (w opt ) + an (w opt ) 



a (w op t) 



Aa (w opt ) = ay (w opt ) - a_i_ (w op t) 



(24) 



(25) 



(26) 



(27) 



Note that all of the components of the polarization vector e in Eq. ( |25j ) are given in the space-fixed 
frame. For CH3F, the static values of the polarizability tensor have been determined experimentally to 
be a (0) = 2.540 A 3 02) and Aa (0) /d (0) = 0.33 E3]. The average polarizability d is smaller than the 
~ 150 A 3 values typical of the alkali metal dimers ll44l l45l . As the electronic structure of CH3F lies in 
the deep ultraviolet [46 1, optical trapping at the common frequencies of 1064nm or 532nm appears to be 
reasonable [47 1. 

For fields whose polarization is parallel to a spherical unit vector, the factors T±2 and T±i all vanish 

identically, and so M remains a good quantum number. However, even in this case states with different 

(2) 

\M\ experience different depths of the optical potential due to the tensor shift depending on Cq (w). As 
the typical depth of an optical potential (< 10-100 kHz) is much less than the splittings between rotational 



levels(^GHz), we can consider only the diagonal terms of Eq. (24 1, where the expected tensor shift is 
proportional to 

3K 2 ] [J(J+1) -3M 2 ] 



(JKM\C^ 2) (lu)\JKM) = 



[J (J 



1) 



J(J+1) (2J-1) (2J + 3) 



(28) 



We expect that many of the techniques which have been devised to cancel the tensor shift in linear rigid 
rotors [28] can also be adapted to the case of symmetric top molecules. 

We now turn to the rotational Zeeman effect, which is described by the Hamiltonian 



H] 



B;rot 



(29) 



where /ijv is the nuclear magneton. The tensor g describes the anisotropy of the coupling between rotation 
and the external magnetic field. As was also the case for the dynamical polarizability above, the tensor g 
has the same symmetry as the inertia tensor, and may be fully described by two elements gii and g±. The 
parameters for CH3F, gii = 0.265 and g± = —0.062, have been determined experimentally [98 1. Taking a 
similar approach as for the AC Stark shift above, we find that the Hamiltonian Eq. (|29]> may be written as 



B;rot 



-Mat 



gB ■ J + Ag £ ( 6 



'7p 



p=-2 



where 



7±2 



Bifi Jzpi 



7±i 



a/3 



70 = - [2B J + B!J_! + B-xJi 



(30) 



(31) 
(32) 
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are linear combinations of rotation operators in the space-fixed frame and B p = B • e p are the components 
of the magnetic field also in the space-fixed frame. Also, we have defined the tensor invariants 

2.9_l + 3|| 

9 = ~ — 3 , A.g = - g± . (33) 

As the nuclear magneton ^ijv ~ 762 Hz/Gauss, the rotational Zeeman effect will not appreciably mix 
rotational levels for any reasonable value of the magnetic field, and so we can consider only the matrix 
elements diagonal in J. Assuming that the field is along the ez direction, B = Bez, we find the matrix 
elements 

K 2 



(JK'M'\H B , rot \JKM) = -fi N BM 



g± + A? 



J(J + l) 



tmm'$kk' ■ (34) 



3 Relationship between symmetric top molecules and magnetic dipoles 

In this section, we first review the single-particle physics of a magnetic atom by showing the explicit forms 
of magnetic dipole moments and the Zeeman effect. Then, we show that there is a rigorous correspondence 
between a magnetic atom and a symmetric top molecule which we present as the mapping in Table [3] 
Finally, we extend our correspondence to include the dipole-dipole interaction, and show that dipole- 
dipole interactions are significantly larger for the simulated magnetic dipole than for naturally occurring 
magnetic dipoles. For example, the symmetric top molecule CH3F can simulate the magnetic atom 52 Cr 
with a factor of 620 enhancement in the dipole-dipole interaction energy. The correspondence shown in 
this section can be applied to a large class of quantum magnets, such as magnetic atoms, long-range spin 
models, and molecular magnets, as will be discussed in Sec. [4] 

3.1 Magnetic dipoles and the Zeeman effect 

Let us consider, as a model of a magnetic atom, e.g., 52 Cr, 168 Er, or 164 Dy, a point particle with an angular 
momentum F and a magnetic dipole operator fi which is proportional to this angular momentum 

A = 9f^bF, (35) 

where gp is the effective g-factor. The appropriate basis for describing this object is \FMp), where 
F (F + 1) is the eigenvalue of F 2 and Mp is the eigenvalue of the projection of F on the space-fixed 



Z-axis. The details of the connection between the effective description of the magnetic dipole, Eq. (35 1, 
and the microscopic angular momentum structure of a highly magnetic atom are presented in Appendix B. 
The magnetic dipole moment couples to an external magnetic field B via the Zeeman Hamiltonian, 

#Zoeman = ~ A ' B . (36) 

The matrix elements of the magnetic dipole operator can be evaluated using the Wigner-Eckart theo- 
rem 1 30 1 as 

(F M' F \jl p \F M F ) = (-lf- M '» ( _ F M ^ * F f ) 5fMb <F||F||F>, (37) 

where the reduced matrix element is 

(F\\F\\F) = ^F (F + 1) (2F + 1) . (38) 

In particular, this gives the first-order Zeeman shifts for a magnetic field along the space-fixed Z-direction 
ez, i.e., B = Bez as 

(FM F \H ZccmaD \FMp) = -gpMpiiBB. (39) 

Hence, insofar as the Zeeman effect is not of the same order of magnitude as the distance between adjacent 
F-manifolds, we find that the eigenstates of a magnetic dipolar particle in a magnetic field are \FMp) 



with energies given by Eq. ( 39 1 and dipole matrix elements given by Eq. ( 37 1 
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3.2 Mapping between magnetic dipoles and symmetric top molecules 

We are now in a position to demonstrate a mapping between symmetric top molecules and magnetic atoms 
which allows us to perform quantum simulation of a gas of the latter with the former. We begin by col- 
lecting the main results from Sees. [5] and 3.1 To first order in cIEdc/Bq, the eigenstates of a symmetric 
top molecule in an electric field Edc = ^DC e z are the states | J KM) with energies and dipole moments 
given by 



(H Iot +H Stalk )\JKM) = 



B J(J+l) + (A Q -B )K 2 - 



dK 



J(J + 1) 



ME] 



DC 



(JKM'\d p \JKM) = (-1) 



J-M' 



J 1 J 

-M' p M 



dK, 



' 2J+1 
J ( J + 1) 



\JKM) , (40) 
(41) 



For a fixed (J, i-Q-manifold, the rotational energy Bo J (J + 1) + (Aq — Bq) K 2 is common to all states 
and so may be taken as the zero of energy. For a magnetic dipole with angular momentum F in a magnetic 
field B = Bez, the eigenstates are \FM F ) with energies and magnetic dipole moments given by 



Hzeeman\ FM F ) 
(FM'p\flp\FM F ) 



g F fj, B M F B\FM F ) , 

F-M' I F \ 



(-1) 



F 



—M'p 



p M f 



(42) 

9F^B^F (F + 1) (2F + 1) . (43) 



Comparing the electric dipole matrix elements Eq. (41 1 with the magnetic dipole matrix elements Eq. (43 1, 
we see that the symmetric top molecule has an electric dipole moment which behaves exactly as a magnetic 
dipole moment provided that we make the identifications 



J ->■ F. 



M 



Mt, 



dK 



F(F+1) 



9f^b 



(44) 



We use the notation "— >" to denote the mapping, as the two sides clearly do not have matching units. 
What we mean in the mapping is that the two systems obey the same algebra provided that we identify the 



coefficients according to Eq. ( 44 1. Now, comparing the energy of the symmetric top molecule Eq. ( 40 1 in a 



fixed (J, if )-manifold and the energy of the magnetic dipole Eq. (42 1 while applying the mapping Eq. (44 1 
we find that the energies are the same provided we make the identification 



Edc 



B. 



(45) 



Thus, a symmetric top molecule in a fixed (J, if )-manifold and a magnetic dipole have the same single- 
particle energies and dipole moments, when the dipole moments and fields are mapped according to 
Eqs. ([44} and ( [45) . 

To complete the mapping, we will define the relationship between the interaction of two symmetric top 
molecules in the same (J, If )-manifold and the interaction of two magnetic dipoles. For this purpose, it 
is useful to recast the dipole-dipole interaction Eq. ([T| as a contraction of two rank-two spherical tensors 
as EE] 



DDI;X> — 



V6Cy 

4 7rr 3 



£ (-If ci>,, 



p=-2 



(2) 
P 



(46) 



where 6 and <fi are the polar and azimuthal angles between particle 1 and particle 2. We have defined the 
irreducible tensor product of two vector operators A and B as 

1 



A®B 



(2) 
V 



1 



1 



m p — m 



A R 



(47) 
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Table 3 Correspondence between symmetric top molecules in a fixed (J, if)-manifold subject to a DC electric field 
and magnetic dipoles in a magnetic field. Note that the correspondence in Table[3]for the special case of a linear rigid 
rotor, i.e., K = 0, gives a vanishing effective dipole moment. That is to say, there is no such correspondence between 
magnetic dipoles and linear rigid rotors. 





Magnetic dipole 


Symmetric top molecule 


Angular momentum 


F 


J 


Effective dipole moment 


gF^B 


dK/J(J+l) 


external field 


B 


Euc 


Dipole-dipole coupling coefficient 


Mo 


1/eo 



Let us consider taking matrix elements of Eq. (46 1 for two magnetic dipoles. Using the basis \FM\M?) = 
\FMx) ® \FM 2 ), we have 

(FM[ M' 2 1 Hddi ; (j, | FM\ M 2 ) (48) 



30mo ( _ 1)2 f-m;-a^ f{f + 1) {2F + 1} g 2 plJ 2 B 



V (Q <h\ V ( 1 1 2 \( F 1 F \( F 1 F 
-f ^' ^ 2^ I m p - m -j) II -Mj m Mi J i -Mg p - m 



x 

p— — 2 m— — 1 



M 2 

where we have used Eq. ( |43) to evaluate the dipole matrix elements. Similarly, for two symmetric top 
molecules in the same ( J, if)-manifold we can evaluate the matrix elements in the basis \ JKM\M2) = 
I JKM 1 ) <g) I JKM 2 ), finding 

(JKM[M^\H DDl . d \JKM 1 M 2 ) (49) 
4ne r 3 y ' J ( J + 1) 



x E 

p— — 2 rn — — 1 



w ^ v f 1 1 2 H J 1 J H J 1 J 

^ l to p - m -p J \ -M{ to Mi M -M^ p - m M 2 



where we have used Eq. ( |4T] > to evaluate the dipole matrix elements. Now, using the mapping Eq. ( |44| , we 
find the mapping between the dipole-dipole matrix elements 

(JKM{M 2 \H nDI . d \JKM 1 M 2 ) -> c 2 <FM{M^|ff DD i ;M | J FM 1 M 2 ) , (50) 



where c 2 = l/(eoMo) is the speed of light in vacuum. The relations Eqs. (44i, (45 1, and (50i establish a 
complete correspondence between the one- and two-body matrix elements of symmetric top molecules in 
a fixed ( J, if )-manifold subject to a DC electric field and those of magnetic dipoles in a magnetic field, 
and hence between the many-body Hamiltonians of dilute gases of symmetric top molecules and magnetic 
dipoles. This correspondence is summarized in Table [3] and may be seen visually in the special case of 
J = 1 in Fig. [3|c)-(f). Appendix C gives a complementary discussion in terms of the second quantized 
Hamiltonians. 

To stress the usefulness of the mapping in Table[3] we will discuss dimensional quantities. The mapping 
in Table [5] then tells us that the ratio of dipole-dipole interaction energies of a magnetic dipole of spin F 
simulated with a symmetric top molecule and a true magnetic dipole of the same spin with <?-factor gp is 

Z?d j2 ts2 

^dd _ fl D^ x 1 2 >< IP 4 C51i 

where do denotes the permanent electric dipole moment of the symmetric top molecule in units of Debye. 
For example, if we were to perform quantum simulation of 52 Cr [5] with CH3F in the J = K = 3 manifold, 
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then the relevant parameters are gp = 2, F = 3, and d = 1.85, giving an enhancement of w 620 in the 
dipole-dipole interaction energy. Even for highly magnetic atoms, such as 168 Er in the F = 6 state [9 | and 
164 Dy in the F = 7 state ifTOl . we see enhancements of 600 and 310, respectively. Thus, a quantum 
simulation of magnetic dipoles with symmetric top molecules allows us to access the physics of magnetic 
dipoles with much larger dipole-dipole interaction energies than are available with even the most magnetic 
atoms. 

A few other notes regarding the mapping in Table [3] are in order. First, a very powerful feature of 
the mapping is that one can choose the effective spin of the simulated magnetic dipole by choosing the 
rotational state of the symmetric top molecule. Hence, while each atomic species has a fixed maximum 
angular momentum, one can choose in principle any effective spin degree of freedom using a symmetric top 
molecule. The only restriction on the effective spin is that it be an integer, as the rotational quantum number 
J is always an integer. This does not imply that the simulated magnetic dipoles are bosons, however, as 
the effective spin of the simulated magnetic dipole is not associated with its statistics. Indeed, symmetric 
top molecules have hyperfine structure, as is discussed in Sec. [5] and the hyperfine degrees of freedom 
will determine the statistics of the simulated magnetic dipole. Of the molecules in Table [2] 13 CH 3 F and 
12 CHaCN are fermionic while the others are bosonic, but this list is by no means exhaustive. The fact that 
the effective spin and statistics of a symmetric top molecule do not have to be related opens the avenue of 
studying magnetic dipoles with integer angular momentum but fermionic statistics, something which does 
not occur in nature with elemental magnetic dipoles. Some applications of our mapping for many-body 
physics will be discussed further in Sec.|4] 

4 Applications of symmetric top molecules for many-body physics 

In this section, we discuss four applications of our correspondence between symmetric top molecules 
and magnetic dipoles for many-body physics, each of which can be called a quantum simulator. The 
first example is to simulate a gas of polarized dipoles. Symmetric top molecules can realize such a gas 
without demanding a large electric field as in 1 E molecules but with hundreds of times larger dipole- 
dipole interaction energies than those of magnetic atoms. Second, by loading a deep optical lattice with 
symmetric top molecules, we can obtain lattice spin models with long-range interactions. Any integer 
spin system is available by choosing the corresponding rotational J-manifold of symmetric top molecules. 
Third, simulation of molecular magnets is another example. Here, we suggest the possibility to simulate 
effective spin models for complex crystalline compounds with a gas of symmetric top molecules in an 
optical lattice where the effects of disorder are more readily controllable. Finally, we present an application 
for the multi-component physics in magnetic atoms. In particular, we discuss the Einstein de-Haas effect 
and spin relaxation dynamics where the dipole-dipole interaction plays an important role and can be more 
clearly observed with the use of symmetric top molecules. 

4. 1 Polarized dipoles 

The simplest scenario of many-body physics with dipolar particles is a gas of single component, point 
particles which possesses resonant dipole moments. Several novel physical phenomena occur for such 
a system, including the manifestation of a roton mode in a dipolar Bose-Einstein condensate l48l l49l . 
BCS pairing in a dipolar degenerate Fermi gas [50, 51 1, crystalline phases in the two-dimensional strongly 
coupled limit |52 53] [54), supersolid [55, 56 571 EH an d topologically ordered phases |59| for dipolar 
gases in optical lattices, and vortices in rotating dipolar gases ||60][6T]|62l[63][64], t0 name a few. More 
information about the physics of polarized dipolar gases may be found in a number of recent reviews of 
the subject ll65l l66l l67l . 

In order that a multicomponent particle realizes a single-component polarized dipole, three conditions 
must be met: (i) a single internal state is occupied, (ii) this state should possess a resonant dipole mo- 
ment, and (iii) all other internal states should be separated by energies large compared to the characteristic 
dipole-dipole interaction energy such that state-changing collisions do not occur. For linear rigid rotors, 
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the condition of a single internal state is easily satisfied by populating only a single level, say, the rota- 
tional ground state |000). In the linear rigid rotor case, however, large electric fields must be applied to 
polarize a significant fraction of the dipole moment in the space-fixed frame, see Table [T] When any rotor 
eigenstate has an appreciable resonant dipole moment due to application of a static electric field, all other 
dipole-allowed states lie naturally very far away in energy. For realizing the single-component polarized 
dipole with magnetic dipoles we assume that only the maximally stretched spin state, \FMp = F), is 
populated, as this is the state with the largest dipole moment. Then, applying a magnetic field such that the 
Zeeman splitting between states -Efi c id = 9f^nB 3> -Edd, where £"dd is the characteristic dipole-dipole 
interaction energy, requirement (iii) is met. Often, the magnetic fields required are very small for mag- 
netic dipoles. For example, in Chromium, a field on the order of 1 mG is sufficient |68 1. For symmetric top 
molecules, we assume that only the | J J J) level is populated, as it has the largest dipole moment. Applying 
an electric field such that -Efieid = ^-Edc ^ £dd> we again satisfy requirement (iii). The fields required 
are on the order of V/m, much smaller than the kV/cm fields required to polarize linear rigid rotors. How- 
ever, larger fields may be applied, and will result in larger effective dipole moments due to (usually small) 
Stark effects beyond first order. For J = K = 1, already we obtain half of the resonant dipole moment. 



For higher J = K states, we can obtain a larger fraction of the permanent dipole moment, see Eq. ( 13 1. 
Hence, we can obtain gases of single-component polarized dipoles with symmetric top molecules which 
interact orders of magnitude more strongly than magnetic dipoles and require orders of magnitude smaller 
electric fields than linear rigid rotors. 

4.2 Lattice spin models 

Let us now consider the case in which a gas of symmetric top molecules in a specific (J, K) -manifold has 
been loaded into an optical lattice such that there exists exactly one molecule per lattice site. We note that 
trapping of molecules in an optical lattice is demonstrated technology for 1 S alkali metal dimers; KRb has 
been trapped in an optical lattice [69|, and RbCs is made directly in optical lattices lfT31l70ll . to name two 
examples. In the limit that the lattice is very deep, which can be achieved by increasing the intensity of 
the optical lattice beams, all molecules reside in the lowest band of the lattice and the motional degrees of 
freedom are completely quenched. Hence, the dynamics of the system reduces to that of the dynamics of 
the internal degrees of freedom. For simplicity, we assume that the tensor shift, Eq. ( |28] i, has been canceled, 
see Sec. |2.3| Then, the many-body dynamics of the internal degrees of freedom may be described, up to 
constant terms, by a lattice spin model 

H Spin - yEE £ A?;™SrS?- m -hY,S!, (52) 



with 



2 

p— — 2 m— — 1 



d 2 K 2 dKE DC 

~ W*3 J2 (J + 1) 2' ft -J(J+l)' W 



A\f ee -730 ( 1 1 2 ) /drdr'K(r)| 2 ^ V ^ \ Wj (r')f . ,54, 
\ to p — m —p J J r — r" 



(2) 

Here, a is the lattice spacing, the functions Cp (8, </>) are the unnormalized spherical harmonics defined 



in Eq. (26 1, Wi (r) is a lowest band Wannier function [71 1 centered on lattice site i, and we have defined 



the spin- J operators of site i, Sf, by the matrix elements 

(JKM'\Sf\JKM) = (-l) J - M ' ( J M , J J M \ y/j(J + 1)(2 J + 1). (55) 



The operators Sf mimic the behavior of a local spin of effective magnitude S = J according to Eqs. (37 1 



and (38 1. The model Eq. (52 1 describes elemental spins of a length fixed by the rotation of the symmetric 
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top molecule interacting via dipole-dipole interactions with characteristic coupling energy W and placed 
in an effective magnetic field h. The coefficients A?'™ describe the anisotropy of the process in which the 
net Z-component of the spin changes by p as a function of the distance between sites i and j. Neglecting 
effects due to anisotropic confinement ll72"ll . the spatial dependence of A?'™ is given approximately as 

C_p(9ij, 4>ij)\i — j|~ 3 , where and Vj are the spatial coordinates of lattice sites i and j, respectively, and 
6ij and cj)^ are the polar and azimuthal angles between lattice sites i and j. 



A limiting case of the spin model Eq. ( 52 1 which is of particular interest is when the effective magnetic 
field h is much larger than the spin-spin coupling W. In this case only the p = processes are resonant, 
and so the total magnetization S?) is conserved. Thus, the magnetic field term may be removed from 
the Hamiltonian by a gauge transformation and we are left with a long-range XX Z model 

i-^j i ^ j i ^ y 

where we have switched from the spherical tensor notation for the spin operators to the more conventional 
ladder operator notation used in quantum magnetism 

S\ = ^ = Sf ±i§y = TV2S? 1 . (57) 

For CH3F in a 532nm optical lattice, the coupling constant W becomes on the order of 15 kHz or larger, 
which is an order of magnitude larger than the corresponding coupling constant for the linear rigid rotor 
proposals with KRb and of the same order of magnitude as the coupling constant for LiCs in the same 



proposals 112711 . An important new feature of our XXZ model Eq. (56 1 compared to proposals using linear 
rigid rotors [27, 28, 73, 29, 74]|75]|76l is that any integral spin S can be achieved in principle without the 
fine-tuning of any fields by simply populating a different rotational manifold J. Systems with higher spin S 
should more closely match the predictions of spin- wave theory, and so our results provide the possibility to 
assess the accuracy of approximate theoretical approaches [78]. We expect that the addition of microwave 
fields to couple rotational states with the same K can greatly extend the tunability of spin models such as 
Eq. ( 52 1, by analogy with what has been done for the case of linear rigid rotors [ 27 , 28 771 125] . 



4.3 Molecular magnets 

One particular realization of large effective spins in condensed matter is crystals of molecular magnets. 
Molecular magnets are organic molecules which contain transition metals with strongly exchange-coupled 
spins, and so each molecule behaves as a single large effective spin [79 1 . Prominent examples of molec- 
ular magnets are the manganese complex Mni20i2(CH3COO)i6(H20)4, referred to as Mni2 l80l IBTI . 
the manganese acetate complex Mni20i2(CH 3 COO)i 6 (H 2 0)4-CH3COOH-4H20, referred to as Mni 2ac , 
and the iron cluster [(tacn)eFeg02(OH)i2] 8+ , referred to as Feg. Here tacn is the cyclic organic ligand 
compound triazacyclononane with chemical formula CgH^NHb. All of the listed molecular magnets 
have effective spins S — 10. However, a range of other effective spins are possible, including half-integer 
spins [79|. Both Mn 12 and Fe 8 may be synthesized as large crystals; for example Mn 12 crystallizes in a 
tetragonal structure. Additionally, as exemplified by Mni2 ac , Mni2 and other molecular magnets can often 
be dissolved in a variety of solvents by ligand substitution, and this can lead to a variety of other crystalline 
structures. 

Molecular magnetic crystals have been of great interest due to the possibility of observing resonant, 
coherent tunneling of magnetization. In particular, the hysteresis loop describing the behavior of the mag- 
netization as a magnetic field is decreased and then increased shows step-like features ll80ll . These step-like 
features have been associated with transfer of magnetization between magnetic sublevels made degener- 
ate by the applied magnetic field, either through thermally assisted [ 80 1 or purely quantum mechanical 
means ll83l . In crystals of molecular magnets, the lattice spacing is large enough that dipole-dipole interac- 
tions between neighboring spins are assumed to be negligible. Instead, tunneling of magnetization occurs 
through terms such as the rhombic zero-field splitting term E(S^. — Sy) which change the magnetization of 
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a single molecule. In contrast, using symmetric top molecules to perform quantum simulation of a crystal 
of molecular magnets, magnetization may be exchanged by strong dipole-dipole interactions. This may 
lead to a richer class of magnetic tunneling phenomena than is seen for conventional molecular magnets, 
for example significant deviations from the predictions of Landau-Zener theory during a magnetic field 
quench 

A quantum simulation of crystals of molecular magnets with symmetric top molecules is possible using 
a setup similar to Sec. |4.2| with symmetric top molecules loaded into an optical lattice. One advantage of 
using symmetric top molecules versus actual molecular magnets is the possibility of strong dipole-dipole 
coupling, as mentioned in the last paragraph. We note that the dipole-dipole interaction can also be made 
weak while the effective spin is still large by choosing a small value for K but a large value for J, see 



Eq. ( 13 1. Another advantage, also discussed in Sec. 4.2 is the ability to change the effective spin S of the 



simulated molecular magnet via the rotational state of the symmetric top molecules. Furthermore, altering 
the optical lattice geometry and the filling allows for control over the crystalline geometry and disorder. All 
of these additional controls may be useful in designing and optimizing molecular magnets for applications. 
Finally, using the tensor shift of the optical lattice or microwave fields, it may be possible to simulate an 
axial zero-field splitting DS 2 Z or the rhombic zero-field splitting mentioned above, both of which play a 
significant role in conventional molecular magnets. 

4.4 The Einstein-de Haas effect and spin relaxation 

Let us now consider a multi-component Bose-Einstein condensate (BEC) with dipole-dipole interactions 
[85). The dipole-dipole interaction Eq. ([T]i does not conserve the projection of the spin angular momentum, 
but conserves the projection of the total angular momentum, where the total angular momentum is the sum 
of internal spin and orbital angular momenta. Hence, the process of the spin relaxation [68 2J]|22], in 
which the projection of the total spin is altered via the dipole-dipole interaction, is associated with transfer 
of angular momentum from the spin to the orbital angular momentum. This transfer of angular momentum 
means that an initially polarized gas of dipoles starts to rotate spontaneously, a phenomenon known as the 
Einstein-de Haas effect [ 86 1 , possibly producing non-singular vortices in the multi component BEC. Based 
on mean field theories [87, 88 1, the dynamics of spin relaxation is interpreted as the Larmor precession of 
dipoles around the effective magnetic field produced by inhomogeneous condensates of bosonic dipoles. 
Thus, the dipole-dipole interaction determines a typical time scale for the spin relaxation, which for bosonic 
atoms is on the order of r sr ~ hjicddP) where n is the density and Cdd = Ho(9fHb) 2 /(47t). For a 
52 Cr BEC with n ~ 10 14 cm -3 , r sr is a fraction of a millisecond, while for symmetric top molecules 
with J = K = 3 it can be shortened to the order of microseconds for the same density due to a six 



hundred-fold increase of Cdd with respect to 52 Cr, see Sec. 3.2 Likewise, spin relaxation dynamics for 
symmetric top molecules can be seen on the same timescales as for Chromium with densities a factor 
of six hundred smaller. In addition, in order to see spin relaxation, the characteristic energy difference 
between single-particle states with neighboring spin projections M and Mil should be small compared 
to the characteristic dipole-dipole interaction energy. For 52 Cr, the dipole-dipole interaction is weak, and 
so spin relaxation requires very small magnetic fields of 1 mG or less [89 1. For symmetric top molecules, as 



discussed in Sec. 2.2 the energetic differences between states with neighboring projections M and M ± 1 
can be small compared with the dipole-dipole interaction for a wide range of the applied DC electric field. 
There are many other studies on the physics of multi-component magnetically dipolar atoms, such as spin 
textures [90][9T||2l and dipolar relaxation [19, 20 1; we expect that quantum simulation with symmetric top 
molecules can access these phenomena on shorter timescales or with lower density and even explore new 
regimes with large dipole-dipole interaction energies compared to isotropic contact interaction energies. 

5 The influence of hyperfine structure 

The correspondence between symmetric top molecules and magnetic dipoles has been discussed within the 
context of the rigid rotor approximation (RRA), where the nuclei are structureless objects rigidly connected 

Copyright line will be provided by the publisher 



18 



Wall, Maeda, and Carr: Simulating quantum magnets with symmetric top molecules 



in space. Within this approximation, the states | J KM) and \J, —KM) are degenerate due to the cylin- 
drical symmetry of the molecule, and possess resonant space-fixed dipole moments ( JKM\do\JKM) = 
dKM/[J( J+l)] even in the absence of an electric field. As mentioned in the introduction, the existence of 
an elemental electric dipole moment would violate both parity and time reversal symmetries J99j. While 
neither parity nor time reversal are good symmetries of the universe within the context of the standard 
model, their violations are extraordinarily slight [2 100|. Thus, it is expected that the existence of space- 
fixed dipole moments for symmetric top molecules in the RRA is a defect of the chosen approximate basis 
rather than being due to breaking of near-fundamental symmetries of the universe. Any interaction which 
couples the states | J KM) and \J, —KM) will result in the proper physical eigenstates in the absence of 
fields being instead the linear combinations | J KM) ± \J, —KM), and introduce some energetic splitting 
Ajkm between these states. These linear combinations have no expected dipole moment, and hence will 
not display a linear Stark effect. The transition to the basis in which | J KM) and | J, —KM) are decou- 
pled and hence a linear Stark effect is observed occurs for electric field strengths of order AjxAi/d- In 
this section, we will show that the splittings Ajkm are on the order of a few kHz for the \K\ — 1 levels 
and un-observably small for all higher \K\ levels. Hence, the linear Stark effect holds for essentially any 
non-vanishing field for \K\ > 1, and for fields greater than V/m for \K\ = 1. The fields where the linear 
Stark effect begins to break down are many orders of magnitude larger than this, see Table [2] In this sense, 
accounting for terms which couple together states with ±\K\ does not significantly affect the mapping of 
See. [3] 

For molecules with C^ v symmetry]^ and non-zero nuclear spin such as CH3F, states with \K\ ^ 3N (N 
an integer) are degenerate to all orders in the rotation-vibration Hamiltonian. Thus, breaking of the ±\K \ 
degeneracy dominantly occurs through hyperfine effects for such molecules. In the following subsections 
we will investigate the role of hyperfine structure in the rotational structure of the canonical symmetric 
top molecule CH3F, shown in Fig. [T] A more general discussion of the breaking of the ±\K\ degeneracy 
for non-planar symmetric tops of the form XY3, including symmetric tops with no nuclear spin, may 
be found in Ref. H 1 II - The dominant hyperfine interactions of symmetric top molecules without nuclear 
quadrupole coupling are spin-rotation interactions and spin-spin interactions. We discuss these interactions 
in Sees. 5.2 and 5.3 respectively. Before discussing these interactions, we briefly review the classification 



of molecular levels according to point group and exchange symmetries of constituent nuclei in Sec. 5.1 
for non-specialists. Finally, in Sec. |5.4| we discuss the nuclear Zeeman effect. In addition to discussing 
hyperfine interactions for the purpose of clarifying the degree of breaking of the ±\K\ degeneracy, isolation 
of particular hyperfine states is often an important step in reaching quantum degeneracy l34l . and so we 
compile information about the hyperfine structure of CH3F here to aid future research on the many-body 
physics of symmetric top molecules. 



5.1 Classification of molecular symmetry 

A classification of the symmetry of a molecule using the language of point groups results in the full 
permutation-inversion group of the molecule [102] in which the eigenstates of the molecule transform 
irreducibly under spatial inversions as well as exchanges of identical nuclei. Using group theory, it can be 
shown that the only states allowed by exchange symmetry belong to one of two non-degenerate represen- 
tations which have opposite parity II 1031 11021 . For CH3F, the molecular symmetry group is C^ v , which 
corresponds to a 3-fold rotation axis which exchanges the three Hydrogen nuclei in a clockwise fashion as 
well as reflections through any plane which contains the symmetry axis and one of the Hydrogen nuclei. 
There exist three representations of Cs v , denoted A\, A2, and E. The representation A\ is completely 
symmetric under all group operations, and both A\ and A2 are non-degenerate. The representation E is 
doubly degenerate. In addition to the rotational wavefunctions, the three Hydrogen nuclear spins can cou- 
ple to form either ortho-CH 3 F (7 H = 3/2) or para-CH 3 F (J H = 1 /2), with J H the total nuclear spin of the 



A brief review of molecular point group symmetry is given in Sec. 



5.1 
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H nuclei. A detailed examination of the molecular symmetry group, taking into account indistinguisha- 
bility of the Hydrogen nuclei, demonstrates that not all combinations of nuclear and rovibrational wave 
functions give internal wave functions with the appropriate symmetry. In particular, when K = 3N, N a 
positive integer, the rovibrational wave function has symmetry A\ © A2, and so the nuclear spin function 
must belong to the completely symmetric representation A±, ortho-CH 3 F, in order that the internal wave 
function has the correct symmetry. In contrast, when K = 3N ± 1 the rovibrational wave function has E 
symmetry, and so the nuclear spin wave function must also have E symmetry, para-CH 3 F, in order that the 
complete wave function transforms properly. The actual rovibrational and spin wavefunctions correspond- 
ing to these representations have been given in the literature Il37lll02lll04lll05l . and the methods of group 
theory can be applied to simplify evaluation of matrix elements, just as for the case of SU(2) 1 106 1. 

5.2 Spin-rotation interactions 
Nuclear spin-rotation interactions may be written as 

#s-r = ~J2J2 iX ™ < 58 > 

X i 

where X is a label for nuclei, e.g., X — {H,C,F} for CH3F, and i runs over the different members of 

-X, 



that nucleus, e.g., takes H(i), Hp), H( 3 ) for CH 3 F. The tensor element c a p expresses the magnetic 
field in the a th direction at nucleus Xu\ due to rotation of the molecule about direction (3. For a nucleus 
X which lies on the symmetry axis, the principal axes of the tensor c x are parallel to those of the inertia 
tensor, and c x also has a doubly degenerate eigenvalue. The matrix elements can be computed by analogy 



with the rotational Zeeman effect, Eq. d30b, with the operators Ip w taking the place of magnetic field 



Lp 

projections B p . Hence, as with the rotational Zeeman effect, the spin-rotation interaction for nuclei on the 
symmetry axis will not mix states with different values of K. For CH 3 F, the values = 4.0kHz and 
c| = —51.1kHz have been determined experimentally B107I . and = —6.03kHz and cj? = —18.33kHz 
have been computed for 13 C ( 12 C has zero nuclear spin) 111081 . 

Let us now consider nuclei which do not lie on the symmetry axis, such as the H atoms in CH3F. If 
we define the coordinate system such that the molecule-fixed x-axis points from the symmetry axis to the 
Hydrogen nucleus labeled 1, then reflection symmetry about the plane passing through the symmetry axis 
and this Hydrogen nucleus dictates that the spin-rotation tensor c 11 ' 1 ' takes the general form 

c* y j . (59) 

in the molecule-fixed frame. The tensor c W is not a priori symmetric when computed in a specific 
reference frame, and so the spin-rotation Hamiltonian of the first Hydrogen nucleus must be made mani- 
festly Hermitian 11051 as ff s _ r; H (1) — > \ -ffs-nH (1) + ^l- r H (1) • This amounts to replacing c^ 2 with 

{d^ z + c^ x J /2 and c*} x = (c^ x + c^ 2 *^ /2. Because of the three-fold rotational symmetry of the 
molecule, we have that c H < 2 > = R 2 1 c u wR 2 and c H < 3 ) = R^^^R^, where 



R 2 - R3 1 = I -V3/2 -1/2 I (60) 


is the matrix which rotates the molecule by 27r/3 about the symmetry axis. The values of the Hydrogen 
nuclear spin-rotation tensor have been calculated to be c£ x — —1.25, c^ y — 2.59, c^ z = 15.92, c^ 2 = 5.77, 
= 1.41 for 12 CH 3 F [ 108], with deviations of a few percent for 13 CH 3 F 11081 . As opposed to the case of 
spin-rotation interactions for nuclei on the symmetry axis, the transformation of the off-axis spin-rotation 
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interaction to the space-fixed frame will not result in a summation of spherical harmonics. Rather, the 
summation will involve rotation matrices ^AMAifO^O w ^ both AM and AK nonzero. In particular, this 
implies that such an interaction can mix levels whose values of K differ by at most two, as the interaction 
involves at most rank-two tensors. A physical way to interpret this result is that the Hydrogen nuclear 
spin-rotation interaction breaks the cylindrical symmetry of the molecular Hamiltonian, as the Hydrogen 
nuclear spins have only discrete rotational symmetry about the symmetry axis. The only pair of projections 
K which are degenerate and connected by AK — ±2 are the \K\ = 1 states. Here, the Hydrogen spin- 
rotation interaction splits the symmetric and antisymmetric linear combinations of | JIM) and | J, — 1M) 
by a few kHz, and the linear behavior of the Stark effect does not set in until the field strength becomes 
comparable to this splitting. For states with \K\ > 1, the spin-rotation coupling only acts in second 
and higher order. The resulting splittings between the symmetric and antisymmetric linear combinations 
of | J KM) and \J, —KM) are of the order of 10 _4 Hz or less, and further decrease exponentially with 
increasing \K\ [ 101 1. 



5.3 Nuclear spin-nuclear spin interactions 

Neglecting any electron-mediated interactions between nuclear spins 111 1 1||4T1|109II1 10 1. the nuclear spin- 
nuclear spin Hamiltonian contains only dipole-dipole interactions between all non-zero nuclear spins in 
the molecule. It may be written as 



- V6 /j, H 2 N \ - g x gx 

it. . — 



<x w ,x>> 



(2) 



1(2) 



(61) 



where the notation < X^^Xi^j > denotes a non-ordered pair of X^ and XL,-., that is, each pair of 
nuclei is counted exactly once, R X(i) x' { is the distance between the nuclei X^ and X'^y and ex (i) x' { , 
is a unit vector pointing from nucleus X^ to nucleus XL,^, The nuclear (/-factors gx are well-known and 
tabulated 11 1211 . For convenience, we note here that the g-factors relevant for CH 3 F are gn = 5.585694, 
g F = 5.257736, and gc = 1.404822 for 13 C. Spin-spin interactions between the Hydrogen nuclei only 
contribute when the ortho configuration of nuclear spin is involved, and hence are only present when 
K = 3N with a positive integer N. On the other hand, the Fluorine-Hydrogen spin-spin interactions 
contribute for any configuration of the Hydrogen nuclear spin. Similarly to the Hydrogen nuclear spin- 
rotation interaction, the Fluorine-Hydrogen spin-spin interaction involves rotation matrices ^AMAif ( w ) 
when expressed in the space-fixed frame, and so it mixes states with K values differing by at most AK = 
±2. As discussed in the above, this interaction causes splittings of a few kHz for the \K\ = 1 levels II 1 1 OH . 
but does not have observable consequences for higher values of K. Hence, we conclude that the ±K 
degeneracy is broken slightly for the \K\ = 1 levels, but can be considered to hold for all other levels 
\K\ > 1. In turn, this means that the linear Stark effect which forms the basis of our correspondence 
between symmetric top molecules and magnetic dipoles also holds for these states even in the presence of 
hyperfine structure. 

5.4 The nuclear Zeeman effect 
The hyperfine interaction we lastly consider is the nuclear Zeeman effect, with Hamiltonian 

ffz ; nuc = -MJV E S 9xI X ^ ■ B . (62) 

X i 



The g-factors were listed following Eq. (61 1 above. For a magnetic field aligned along ez, the matrix 
elements are completely diagonal in the space-fixed frame: 

{JKMI n M n M c M F \H z , lmc | JKMI H M H M C M F ) = - f i N (g H M H + g c M c + g F M F ) B , (63) 

where Mh is the sum of nuclear spin projections of Hydrogens. This expression holds for all representa- 
tions of the Hydrogen nuclear wavefunction. 
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6 Conclusions 

We have shown that the many-body physics of symmetric top molecules in precessing states maps onto 
the many-body physics of magnetic dipoles with effective spin set by the rotational quantum number of 
the molecule. The strength of dipole-dipole interactions in the simulated magnetic dipoles is typically two 
orders of magnitude larger than for naturally occurring magnetic dipoles, which opens up new regimes 
of dipole-dipole interacting systems. We illustrated several manifestations of novel many-body physics 
and quantum simulation hence possible with symmetric top molecules. In one example, symmetric top 
molecules can generate a gas of single-component polarized dipoles with either bosonic or fermionic statis- 
tics. The preparation of such a gas requires much smaller static electric fields than a polarized gas of linear 
rigid rotors, and will feature larger dipole-dipole interactions than a polarized gas of magnetic dipoles. 
Furthermore, we showed that symmetric top molecules in a deep optical lattice may be used to engineer 
quantum simulations of lattice spin models with anisotropic, long-range, and magnetization-changing in- 
teractions as well as arbitrarily large spin. Trapping of symmetric top molecules in an optical lattice also 
allows for quantum simulation of molecular magnets, where symmetric top molecules offer the advantage 
over condensed matter realizations of tuning the geometry, disorder, strength of interactions, and effective 
spin. Finally, we discussed quantum simulation of spontaneous demagnetization and the Einstein-de Haas 
effect on much shorter timescales or lower densities than with naturally occurring magnetic dipoles. The 
existence of precessing states is a feature of symmetric top molecules which is not shared by linear rigid 
rotors, such as the 1 S ground states of the alkali dimers. In addition to the correspondence between sym- 
metric top molecules and magnetic dipoles, we also discussed some technologically important aspects of 
symmetric top molecules, such as opto-electrical cooling to reach quantum degeneracy, optical trapping, 
and the effects of hyperfine structure. 

Our work has just begun the study of many-body physics with symmetric top molecules. In particular, 
we have only studied the cases of static or far-off-resonant fields. We expect that the use of microwave 
dressing will greatly enhance the tunability of effective many-body models with symmetric top molecules, 
just as has been the case for linear rigid rotors [27 28 1. In addition, threshold scattering of symmetric 
top molecules at ultralow temperatures is an unexplored area and beyond the scope of the present paper. 
Studies of ultracold scattering will provide insight into the relative role of short-range contact interactions 
to the long-range dipole-dipole interactions considered in this paper. 
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7 Appendix A: Matrix elements for symmetric top molecules 

In this section we compile useful matrix elements for symmetric top molecules. As the wave functions of 
symmetric top molecules are proportional to Wigner D-matrices, see Eq. |6]), most matrix elements can be 
computed using the theorem on integration over products of three rotation matrices [30 1 

f du@ 3 \ (u) (u) @ h , (u) = Stt 2 ( j \ 32 j3 , ) ( 31 32 33 ) . (64) 

In particular, using the facts that ® 3 m0 (w) = (oj) and &„* m , = (-l) m ~ m ' # m _ m „ we find for the 
expectation of any unnormalized spherical harmonic 

(j'K'M'\Cl e) (w) \JKM) 

= (-D M '-V(2J' + 1)(2J + 1)( _ J M , I J M )[ f K I J K ) 5 K , K . (65) 
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As for the expectation of the dipole operator along a space-fixed direction e p , we note that the permanent 
dipole moment of a molecule must lie along the symmetry axis. We define the symmetry axis to be the 
molecule-fixed z-axis e z = e q= o, which yields d = de q= Q. We now rotate from the molecule-fixed frame 
to the space-fixed frame as BTI 

1 

i = E 9%d 9 = dcP. (66) 

q=-l 

Hence, the expected dipole moment along space-fixed direction e p is given by Eq. ( 65 1 as 

(J'K'M'\d p \JKM) 
- (-l) M '-V^< + l)(2 J + l)(_ J M , I J M ){! K J J K ) 



>K'K 



(67) 



8 Appendix B: Microscopic angular momenta and effective angular mo- 
mentum of a magnetic dipole 

In this section, we discuss the relationship between the algebraic structure of the intrinsic angular momenta 
of an atom and the effective magnetic moment associated with the total angular momentum of the atom. 
Let us consider an atom with orbital angular momentum L, spin angular momentum S of electrons, and 
nuclear spin angular momentum I of nuclei. The magnetic dipole operator of the atom is 

fi = \x B (git + g e Sj +fi N g N I. (68) 

Now, we consider the fully coupled representation in which we first couple L and S to form J, J = L + S, 
and then couple J and I to form the total angular momentum F, F = J + I. The associated basis is 
denoted by {\((L, S)J, I)FM F )}. Using standard angular momentum decoupling techniques [30|, we 
can compute the projections of the intrinsic angular momenta in this basis as 

(((£', S')J', I')F'M' F , | L p |((L, S)J, I)FM F ) 

= 5 LL ,5 ss -5ii> (F'M' F ,\F P \FM F ) 

[F(F + 1) + J(J + 1) - /(/ + 1)] [J(J + 1) + L(L + 1) - S{S + 1)] 



2F(F+1) 2J(J+1) 
(((L', S')J', I')F'M' F , | S p \((L, S)J, FjFM F ) 
= Slu&ss'Siv (F'M' f , I F p \FM F ) 

[F(F + 1) + J (J + 1) - 1(1 + 1)] [J (J + 1) + S(S + 1) - LjL + 1)] 
X 2F(F+1) 2J(J+1) 

(((L', S')J', I')F'M' F ,\i p \((L, S)J, I)FM F ) 



(69) 



(70) 



= o LL ,d ss >5ii> (F Mp,\F p \FM F ) 2F(F + l) ' ^ 

Importantly, we see that all three projections are proportional to the projection of the total angular momen- 
tum F as a result of the fact that all intrinsic angular momenta are spherical tensors. Hence, we obtain for 
the projection of the magnetic dipole operator, Eq. (|68]l, 



(((L\S')J\I')F'M'p,\fL p \((L,S)J,I)FM F ) = g Fl i B (F ' M' F ,\F p \FMp) 5 w 5 SS ^ir , (72) 
where the effective g-factor, or generalized Lande g-factor, is given by 

g F = C(F,J,I)[ gi C(J,L,S)+g e C(J,S,L)} + ^g N C(F,I,J), (73) 

Ms 
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with a function 

d*,M) = «(« + '> + W + ' ) -^ + 1) . (74) 

2a(a + 1) 

Since hn / — rn e /m p ~ 5 x 10~ 4 , the last term in the braces may be neglected when considering the 
overall magnitude of the magnetic dipole. However, it is important to note that the nuclear spin implicitly 
contributes to the magnetic dipole moment via the total angular momentum. Also, substituting gi = 1 and 
g e ~ 2, we have 

(((L',S')J',I')F'M F ,\(i p \((L,S)J,I)FM F ) ~ ~g Fl x B (F'M' F , \ F p \FM F ) 5 LV 5 SS ,8 IV , (75) 
with the approximate value for g F , 

g F = C(F,J,I)[1 + C(J,S,L)] . (76) 
When we are focusing on a subspace with fixed values of L, S, J, I, and F, which is denoted simply 



by J£({\FMf)}m =_f)> the effective g-factor Eq. (73 i becomes constant, and thus the magnetic dipole 



moment fi is proportional to to the total angular momentum F as a linear operator in this subspace, 

A = 9FVBF in ^({\FM f )} f Mf= _ f ). (77) 
Note that such a subspace spans an effective Hilbert space in several important cases and that it is useful 



to employ the equality (77 1. The first case is rather trivial where only a single angular momentum L, S, or 



I is present, as then Eq. (77 i becomes true in a whole Hilbert space. This is the case for Chromium, which 
contains only electronic spin in its bosonic isotopes [5 |. In addition, the effective theory holds whenever 
/ is zero and there is large separation between F-manifolds compared to the Zeeman splittings within 
an F-manifold. Such a splitting comes from, e. g., spin-orbit coupling. All of the bosonic isotopes of 
the most highly magnetic atoms Dy and Er have zero nuclear spin due to the fact that they have an even 
number of both protons and neutrons 11921 . a trend which holds for a large majority of the highly magnetic 
atoms. The splittings between F-manifolds are many orders of magnitude larger than any possible inter- 
manifold Zeeman splittings in these atoms 11931 [TOl [5), and so the effective therory works. Finally, the 
effective description holds whenever all three angular momenta are present, but the spacings between 
hyperfine manifolds are large. This is the case for the fermionic isotopes of Dy and Er, as they all have 
large nuclear spins and, accordingly, large nuclear quadrupole interactions on the order of hundreds of 
MHz [94, 1 1 , 93 95 1 . The hyperfine spin-spin interactions in the bosonic isotope 169 Tm [96 1 are similarly 
large, despite the fact that there is no nuclear quadrupole [97|. The description of the magnetic dipole in 
terms of the total angular momentum hence holds up to magnetic fields which are much larger than the 



fields required to reach the regime of polarized dipoles, see Sec. 4.1 Thus, in discussion of many -body 
physics, there is no loss of generality in assuming that the microscopic structure of the magnetic dipole 
moment of a highly magnetic atom Eq. ( |68| l may be replaced by the effective structure Eq. ( |77| . 

9 Appendix C: The mapping from symmetric top molecules to magnetic 
dipoles in second quantization 

In this section, we demonstrate that there is a correspondence between the second-quantized many-body 
Hamiltonians for symmetric top molecules and those for magnetic dipoles. We begin with the second 
quantized Hamiltonian for spin-F magnetic dipoles Q], 



H = /dri4 M >i)#^W P (ri) 

M F = -F'' 

« , "V ~\S ' " 



(78) 



i,j=X,Y,Z 
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where the colons denote normal ordering and ip F M F ( r ) creates a magnetic dipole at position r with its 

The field operator ip FM (r) satisfies commutation or anti- 



3.1 



internal state \FM F ), discussed in Sec. 
commutation relations depending on the statistics of the magnetic dipoles. Here, we have defined the 
single-particle Hamiltonian of the magnetic dipoles , which is composed of the kinetic term and the 



Zeeman term for B = Bz^z, Eq. ( 39 1, 

H[ m) = -|- - g Fl i B M F B z , (79) 
2m 

and the function Qij, the kernel of the dipole-dipole interaction, as Qij{v) — [Sij — 3(e r )i(e r )j]/r 3 . We 
have also introduced local spin operators Mf (r), 

F 

Mf(r) = Yl ^fm' p (*")(/») FM' P ,fMf iPfMf (r) ) (80) 

M' F ,M F = -F 

with the matrix elements of the magnetic dipole moment operator (/$) fm' f .fm f (i — X, Y, Z) defined by 

{fi)FM>FM F = (FM' F ,\f M \FM F ) . (81) 



The explicit form of Eq. (ISTji is given by Eq. (43 i in Sec. 3.2 



On the other hand, the second quantized Hamiltonian for symmetric top molecules is given by 

H = / dr l ^JKM(*M e) TpJKM{Tl) (82) 



J KM ' 

+ \Uv 1 dv 2 ■■A(r,) Qiji 2 r ~ r2) t> j (r 2 ):, (83) 

J i,j=X,Y,Z 

where ^j KM (r) creates a symmetric top molecule at position r with its internal state | J KM), discussed 
in Sec. |2.1| and satisfies commutation or anti-commutation relations depending on the statistics of the 



molecules. Here, denotes the single-particle Hamiltonian of symmetric top molecules, composed of 



the kinetic term and the linear Stark term for E = i?DC e z, Eq. ( 1 1 



H[ e) = -—- ^ , ME DC . (84) 
1 2m J(J+1) DC V ; 



Similar to the local magnetization operators, Eq. ( 80 1, local dipole moment operators T>i (r) are defined by 



>K'M'( r )(di)j'K'M',JKMlpJKM(r) , (85) 

J'K'M'JKM 

with the matrix elements of the electric dipole moment operator for symmetric top molecules (di)j'K'M',jKM 
[i = X,Y,Z), 

{di)j. K >M:,jKM = (J'K'M'\di\JKM) . (86) 
We remark that the summation of J, K, M runs over all non-negative integers for J, and integers such that 



— J < K < J, — J < M < J for K and M. As discussed in Sec. 2.1 for the symmetric top molecules in 
a fixed (J, if) -manifold, we have 

{d P )jKM-.jKM = (-1) J - M ' (_ J M , l p (JK'\d\JK) , (87) 
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with the reduced matrix element given by Eq. (|T5]). Thus, for any given species of magnetic dipoles with 
integral spin-F, we can choose the corresponding (2F + 1) states from rotational states of symmetric top 
molecules with J = F and K such that 

(di)FKM',FKM = {JK 'fl^L {FM' F \h\FM F ) 
g F fj, B {F\\F\\F) 

—^{fi)Fm' F ,Fm F ■ (88) 



g F fi B F(F 



When we concentrate on the fixed (J, i^)-manifold of the Hamiltonian Eq. (82i with J = F and K ^ 0, 
it reduces to 

H = / dr i 4« (ri)Hl e) ip F KM(ri) (89) 

M=-F J 

+ i/dr ld r 2 £ rPf^rx) ^^-^ ^):, 
with local dipole moment operators in the (F, F)-manifold T>f K (r) given by, 

F 

^f K ( r ) = E i'FKM'( V )(di)FKM'.FKM^FKM^) ■ (90) 

M'M=-F 



Comparing Eq. ( |89| ) with Eq. ( [78) , we conclude that the second quantized Hamiltonian for symmetric top 
molecules can simulate the many-body physics of magnetic dipoles. 
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